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A formalism for the relativistic description of hadron decays within the constituent quark model is presented. 

First, hadron amplitudes of the light-cone constituent quark model, in particular the weak transition form 

factors at spacelike momentum transfers, < 0, are represented in the form of the dispersion integrals over 

the hadron mass. Second, the form factors at > are obtained by performing the analytic continuation 

from the region < 0. As a result, the transition form factors both in the scattering and the decay regions 

are expressed through light-cone wave functions of the initial and final hadrons. The technique is applied to 

- - - , the description of the semileptonic decays of pseudoscalar mesons and direct calculation of the transition form 

)C} • factors at (7^ > 0. 
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1 Introduction 



Weak decays of hadrons provide an important source of information on the parameters of the standard model 
of electroweak interactions, the structure of weak currents, and internal structure of hadrons. Hadron decay 
rates involve both the Cabibbo-Kobayashi-Maskawa matrix elements and hadron form factors, therefore the 
extraction of the standard model parameters from the experiments on hadron decays requires reliable information 
on hadron structure. 

^ ' The problem of theory is to describe hadron form factors which involve both perturbative and nonperturba- 

OO . five contributions. Higher order corrections to weak currents are calculable perturbatively and can be predicted 
' to high accuracy. The calculation of hadronic matrix elements of the weak currents inevitably encounters the 
problem of describing the hadron structure and requires the nonperturbative consideration. This gives the main 
' uncertainty to the theoretical predictions for the hadron transition amplitudes, 
l/^ , In the case of the semileptonic K13 decay, the K ^ tt weak transition form factor deviates from unity only at 

On ' the second order in comparatively small S'?7(3)-symmetry breaking and can be calculated to high accuracy 
(— I . that provides the most accurate value of the Vus- For extracting the Vd, Vcs, Vub, and Vcb from the decay rates 
^ 1' and lepton spectra in D12, and B13 decays, a reliable calculation of hadron transition form factors at timelike 
JL , momentum transfers is necessary. 

Q ' In last decade an increasing amount of publications have been devoted both to the perturbative calculation 

, of higher order corrections to weak currents and to the description of weak matrix elements of hadrons. We 
shall concentrate on the latter problem, closely related to the investigation of the nonperturbative aspect of 
hadron structure. 

' Various theoretical approaches have been applied to the calculation of the nonperturbative contribution to 

hadron form factors. The most popular amond them are the quark model [^~|^, QCD sum rules p|-jl5|l, and 
lattice QCD calculations |l^-[|l8j. A critical comparison of these approaches can be found, e.g. in ||20|. 

Systems containing heavy and light quarks are usually considered within the Heavy Quark Effective Theory 
(HQET) an effective theory based on QCD in the limit of infinitely large quark masses. The B D, D* 
decays associated with the heavy-to-heavy (6 c) quark transition are described in terms of a single universal 
Isgur-Wise (IW) function Q which can be estimated with any of the mentioned nonperturbative approaches. 
The 0(l/mg) corrections to this picture can be consistently calculated within the HQET (for a detailed review 
see H). 

For the decays caused by the heavy-to-light quark transitions {D K,K*;D tt,p; and B tt,p) the 
situation turns out to be less definite. The HQET does not work properly in this case, and theory faces at least 
two practical problems, namely: (i) The existing theoretical considerations fail to describe the experimental 
results for the B K*ip decay [Q; (ii) In the absense of experimental information on _B — > tt, p decay modes, 
the uncertainty of the theoretical predictions for relevant form factors is too large to make any definite conclusion 
on their values (see Table |]). 

This stimulates further investigation of hadron transition form factors. 

Our special interest lies in the quark model which reflects at the phenomenological level intuitive ideas on 
hadron structure. Various versions of this model have been used for calculating hadronic matrix elements of 
weak currents. 



*E-mail: melikhov@sgi.npi.msu.su 
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Recently, it has become clear that for a consistent and successful application of quark models to electroweak 
decays, a relativistic treatment of quark spins is necessary pO|] , p3| . However, in the first models by Grinstein, 
Isgur, Scora and Wise (ISGW) §, and Wirbel, Stech, and Bauer (WSB) g quark spins were not treated 
relativistically. A nonrelativistic approach by GISW is based on a successful potential model for meson spectrum 
For the calculation of the electroweak form factors, rescaling of the parameter in the form factor q^- 
depcndence is used. Such an alteration has no strong theoretical deduction. In addition, an extrapolation from 
the truly nonrelativistic region w Qmax^ where the model is rigorously valid, to a highly relativistic point 
(7^ = is performed. The first step to the relativistic treatment of meson decays was done in the WSB approach. 
The quark model calculations are performed only at one point q^ — using the Infinite Momentum Frame. For 
the form factor g^-dependence the authors postulate a monopole behavior determined by the nearest vector 
meson state. Although the model considers the quark motion relativistically, the quark spins are again treated 
in a nonrelativistic manner. The WSB approach, as well as the GISW model and its modifications 
have both the theoretical and experimantal objections, namely: the models do not reproduce the IW scaling 
of the form factors in the heavy quark limit, and fail to describe the data on the widths and polarizations 
in the semileptonic D — {K,K*)lv decays (Table ||). The answer to these difficulties lies in the correct 
relativistic consideration of the spins. The exact solution to this complicated dynamical problem is not known, 
but a simplified self-consistent relativistic treatment of the quark spins can be performed within the light- 
front formalism [^|. A description of electroweak properties of pseudoscalar mesons 0-j9| and transition form 
factors at < [|7| was performed in the framework of this formalism. The only difficulty with this approach, 
is that the applicability of the model is restricted by the condition q^ < 0, while the physical region for hadron 
decays is < < {Mi — M2)^, Afi.2 being the initial and final hadron mass, respectively. So, for obtaining the 
form factors in the physical region and decay widths and lepton distributions, assumptions on the form factor 
behavior were necessary. A procedure to remedy this difficulty is proposed here. 

We present a formalism for the relativistic description of the form factors of hadron decays within the 
constituent quark model. For a direct calculation of the transition form factors at timelike momentum transfers 
we use a dispersion formulation of the light -cone constituent quark model p6| . Namely, the amplitudes of 
hadron interactions considered within the framework of the light -cone formalism are represented as dispersion 
integrals over the hadron mass. After that, the form factors at > are derived by performing the analytic 
continuation from the region < 0. As a result, the form factors in the whole kinematic region q^ < {Mi — 
are expressed through the light-cone wave function of a hadron. The developed formalism is applied to the 
analysis of the electroweak properties of pseudoscalar mesons. 

In the next section we demonstrate the equivalence of the light-cone constituent quark model and the 
dispersion relation approach ||2^, We present all technical details of the description of a pseudoscalar 

meson within the dispersion relation approach (leptonic weak decay, two-photon decay, elastic electromagnetic 
form factor) and show the results to be equal to those of the light-cone quark model. 

Section 3 considers the transition form factors. Firstly, at spacelike momentum transfers the light-cone 
expression is reformulated as a double dispersion integral representation. Secondly, the analytic continuation 
to the region of timelike momentum transfers is performed. Along with the normal Landau singularities, the 
anomalous non-Landau singularities contribute to the transition form factors in this region. 

In Section 4 the electroweak properties of pseudoscalar mesons are considered. The following issues are 
addressed: 

1. The dependence of the axial- vector decay constant fp and the heavy-meson elastic form factor on the heavy 
quark mass toq is analyzed, using a parameterization of the meson wave function based on the heavy quark 
symmetry. The corrections to the leading l/mQ-behavior are estimated to be at the level of 10 20% in the 
region of b- and c-quark masses. 

2. The form factors of pseudoscalar meson decays are calculated. Our results are in agreement with the QCD 
sum rules and the experimental data. The form factors can be approximated to a high accuracy by the dipole 
formula in the physical region, but do not contradict to the vector meson dominance as well. 

3. The correlation between the values of Jo and Jb and the slope of the Isgur- Wise function is studied. The 
parameter is found to be in the range 0.7 < f? < 0.9 for reasonable values of heavy-meson decay constants. 
We also discuss possible reasons of the deviation from unity of the Isgur- Wise function at zero recoil. 

The results are summarized in the Conclusion. The appendix provides relevant technical details of the 
dispersion approach. 
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Figure 1: Meson dispersion loop graph Bps{P^). 

2 Quark structure of pseudoscalar mesons 

An approach to a composite system description based on dispersion relations ^ allows constructing relativistic 
and gauge invariant amplitude of the interaction of a composite system with an external vector field starting 
with low-energy constituent scattering amplitude (see the Appendix A) . Two-particle s-channel interactions are 
consistently taken into account both in the constituent scattering amplitude and the amplitude of interaction 
with an external field. In the case of a bound state, its form factor and structure function are expressed 
through form factor and structure function of mass-shell constituents and the vertex G of constituent-bound 
state transition. This vertex is defined by the two-particle irreducible block of the constituent scattering 
amplitude. On the one hand, the dispersion integral representation turns out to be equivalent to the Bethe- 
Salpeter treatment with a separable kernel of a special form, the vertex G being connected with the amputated 
Bethe-Salpeter wave function of the bound state |29j. On the other hand, this approach is equivalent to the 
light-cone description of a bound state with the special form of spin transformation (the Melosh rotation). The 
vertex G determines the light-cone wave function of the bound state. Because of the relativistic invariance, 
the dispersion integral formulation of the light-cone approach approach does not face the problem of choosing 
appropriate component of the current for calculating the amplitudes of the bound state interaction. 



2.1 The quark— meson vertex 

We discuss the case of a pseudoscalar meson, but the same procedure can be applied to other hadrons as well. 
The pseudoscalar meson P with the mass M is considered to be a bound state of the constituent quark with 
the mass mi and the antiquark with the mass m2. To derive the expressions for the soft amplitudes of the 
meson interactions hke < iii'\P >,< 77I-P0 >, and < P'| J^jP >, we start with the corresponding amplitudes 
of the constituent quark interactions < ^v\QQ >, < "i^\QQ >, and < QQ\Ji_i\QQ > and single out the poles 
corresponding to the meson. The amplitude of the QQ interaction turns out to be the basic quantity describing 
constituent-quark structure of the bound state. Near a bound state with the — 0^, the amplitude is 
dominated by the S'-wave partial amplitude which can be expressed in the two-particle approximation through 
the dispersion loop graph Bps with the vertex 

0°(fci, mi)i75Q''(-fc2, m2) o 



G{P') (1) 



with a a color index, Nc — 3 the number of quark colors, = ml, A;| — and P — ki + k2, P^ = s ^ . 
For on-shell constituents, the expression (Q) is the only independent spinorial structure. 

The dispersion loop graph Figj^, which is connected with the meson vertex normalization, reads 

00 

Bps{P^)= f '^'^V,, Pps{s), Bps{M^)^l (2) 



7r(s-P2)' 

(mi+m2)^ 



with Pps{s) the spectral density of the Feynman loop graph 

Pps{s, mi, 1712) = j dkidk25{kl - m\)5{kl - ml)S{P ~ ki - k2) Sp (^{ki + mi)i"f^{m2 - hjild^ (3) 

_ AV2(s,mf,mi), 



where 



-{s-{mi-m2r)9is-{mi+m2y), 
Sirs 

A(s, m\,m\) = (s + m\ — 'rii\)'^ — Asm\. 
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Figure 2: The series of dispersion graphs for < (5Q|A^(0)|0 >. 

Taking into account constituent-quark rescatterings leads to the renormahzation of G (see the Appendix A), 
and the soft constituent-quark structure of the pion is given by the vertex 

Q°(fci, mi)i75Q°(- fc2, ma) ^ ^^^^ 



where Gy{s) = G{s)/B'{M^), 

Gl{s)pps{s,mi,m2)ds ^ 
7r(s - M2)2 

Once the soft vertex (^ is fixed, we can proceed with calculating meson interaction amplitudes. 

2.2 Weak decay of a pseudoscalar meson 

Let us consider the decay P iiiy. The corresponding amplitude reads 

< P\A';,^0)\0 iP, fp, (6) 

fp is the meson axial-vector decay constant. To obtain the expression for this matrix element we must first 
consider the quantity 

<gg|A^"(o)|o> (7) 

with Q a constituent quark, while the axial current v4^(0) — 9(0)7^759(0) is defined through current quarks. 
Next, we must single out the pole corresponding to the pion. 
The bare matrix element has the structure 

< Q(fci)Q(fc2)|^r(0)|0 >bare= Qih) [7^755° (^'') + Pf.Toh°+{P^) + {h ~ k^) ^-f^h'LiP^)] Qi-k^) (8) 

If current quarks were identical to constituent ones we would have had 

g0(p2)^l, hl{P^)^0, h^_{P^)^Q. (9) 

It is reasonable to assume that at least at — (jyi^ _ 777,2)2 the form factors g and h are not far from these 
values. 

The rescatterings of the constituent quarks lead to the series of the dispersion graphs of Fig.||. The bare 
matrix element enters into a single loop graph whose spectral density is the product of G{s) and the 
corresponding Feynman graph spectral density which reads 



dkidk2S{k( - mi)S{k^ - rni)S{P - fci - fca) (10) 



xSp {b^l59AiP^) + Pt,Toh\{P^) + (fci - ^2)^75/1° (i"')](fci + mi) z75(to2 - ^2)) 
The trace is equal to 

- 4z(fci^m2 + k2^,mi)g\{P^) + AiP^ik^k^ + mim2)h\{P^) + 4i(fci - k2)f,{kik2 + mim,2)h"_(P^) (11) 
So the expression for the loop graph takes the form 

B.-^^P.^N^ J _ M2) Tq^s 2i 

(mi +7712)^ 
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x[(mi +m2)<?!^(p2) - s h\{P^) - {ml - ml)h^_{P^)] 

The amplitude with the quark rescatterings taken into account has the same spinorial structure as the hare 
amphtude 

< Q{k^)Q{k2)\Al-{m >= Q{ki) [i^.1b9a{P^) + P^j5h+{P^) + {h - A;2)^75/»-(P')] ^(-^2) (13) 

with 

The form factor develops a pole at — as Bps{M^) = 1. Near P^ = the pole dominates the 
amplitude 

< QQ\Af,\0 >=< QQ\P > ^ < P\Af,\0 > +regular terms. (14) 
Comparing the pole terms in ([T^ ) and ( |l4| ) and using the relation 

< piQQ >= 

one finds 

<P|A-(0)|0 >=zP^/p (15) 

with 

fp = 4VN,J J^^^; (16) 

x[(mi +m2)50(p2) - s /j"+(p2) - (m? - mDh'LiP^)] 

Assuming that in reality the values of and are not far from the limit (|^), we neglect the terms involving 
and come to the relation 

00 

/p^40^(mi+m.),O(M^) / ^^g.(^) AV2(^X,"^l)(^-(^i-^.)-) (17) 

7 7r(s — Al^) Ibns 2s 

(mi+m2)^ 

2.3 The two-photon decay of the neutral pseudoscalar meson 

We consider the decay of the neutral pseudoscalar meson Pq whose constituent quark structure is described by 
the vertex 

^G.(P^) (18) 
The rate of the decay Pq 27 can be written as 

r = ^a''M^g%^, gp^^ = Gp^^(Af2, 0, 0), (19) 
where the form factor Gp-yj is connected with the amplitude 

< 0| (92)^ar(93)l^' 26„,„3^,^3q2^^gf^Gp^^(M^ g|, g|). (20) 

The electromagnetic current J^'"(0) — g(0)7^g(0) is defined through current quarks, whereas the meson struc- 
ture is described in terms of the constituent quarks. So, for calculating the meson amplitude the constituent 
quark amplitude of the electromagnetic current is necessary. The latter is assumed to have the following struc- 
ture 

<Q{k')\m7^.qmQik)>=Q{k')lMk)fc{q^), q^k'~k. (21) 
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Figure 3: The graph describing the decay P'^ — > 77. 

The constituent charge form factor fdq^) is normahzed such that /c(0) = Cc, the constituent charge. The 
anomalous magnetic moment of the constituent quark is neglected in the expression (H]), but it can be included 
into consideration straightforwardly. 

The single dispersion representation for the form factor Gpjj reads 

/ds G (s ') 
— ^— ^^-^Ap^^(si,S2,S3), S2 = g|,S3=g3 (22) 
7r(si - M^) 

where A p~^^ is determined by the spectral density of the Feynman graph of Fig.|^ 

/ dkidk2dk3SiP -k2- k3)d{k2 - fcl - qs)— 2^^tJ 

X Sp(^ij5im- hha^iim + ki)ja3im + k2)^ = -ea2aal32f33<l2'' la^ ^P-tii^i^ ^2, S3) (23) 

The trace reads 

Sp (^ij^im - fc3)7a2(™ + fci)7Q3(™ + ^2)) = 4mea2a3/32/33g*gf^ (24) 

and we find 



Ap-y^ si,S2,S3 = — —jj- -log — — f-^=====\ (25 

47r Xy^{si,S2,S3) \si- S2- S3- Ai/2(si, S2, S3)\/l - 4to2/s / 



Substituting ( |25| ) into (|2^), one obtains the expression which defines the quantity Gp^-y for off- shell photons. 
For real photons one finds 



47r J 7r(s - M^) s \ 1 _ ^1 _ Am? s / 

2.4 The elastic electromagnetic form factor 

The elastic electromagnetic form factor of a pseudoscalar meson is given by the following matrix element 

< n/kr(0)I^M >= [Pit + Pm)^. F^\q^) (27) 

Pi = P'l = M\ Pm - P'm =q,q^< 0. 

Assuming the following structure for the constituent-quark matrix element of the electromagnetic current 
j-(0) = q-(0)7^g(0), 

2n 



< Q{k[MOh^MO)\Q{h) Q{k[)j^Q{ki) Mq^), (28) 
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Figure 4: The contribution H{q ,mf, 111,2) to the elastic form factor. 

the elastic charge form factor of the meson can be written in the form 

F^'iq^) = h{q^)H{q^,mlml) + f2{q^)H{q\ m^, m?) (29) 

in terms of the form factors H. The quantity H{q^ ,rnl,m2) describes the subprocess when the constituent 1 
interacts with the photon, while the constituent 2 remains spectator. 

The double dispersion representation for the form factor H{q^ ^rril^rn^) (Fig.p reads 

Hiq^mlml) = J X^^) ' ^' gVi, ^1, rn.). (30) 



/ dkidk[dk25{kl - ml)5{k'^ - m\)5{kl - ml)6{P - fci - k2)S{P' - k[ - fca) 



Here A is the double spectral density over and P'^ of the corresponding triangle Feynman graph 
1 

X Sp(^{k[ +mi)7^(fci +TOi)i75(m2 - h)^^^ = 2Pf,{q)A{s',s,q'^\mi,mi,m2) (31) 

with 

P,i<l) = iP-^ P' = P" = iP' - Pf - 

The trace reads 

-Sp {{k[ + mi)7p(fci + mi)75(m2 - ^2)75) 2fci^(s - [mi - 7712)^) + 2fci^(s' - (mi - 7712)^) + 2fc2^g^ (32) 

Multiplying both sides of ( ^l| ) by and using ( p^ ) one obtains at < 

-q^ 

A(s',s,g^|mi,TOi,m2) = ,,^,„, , ^ (s's + (s' + s - q^)m2{mi - 7712) - (?7ii + m2){mi - m2f) (33) 

4A'''^(s , s, g"^) 

^(S - (toi + m2)2)6'(s' - (mi + 7772)2)6* (-^^(.s' + s - ^2 _^ 2(777^ - r772))2 + A(s', S, g2)(g2 _ 4^2)) ^ ^2^0 

with A(s', s, g2) = (g' + s - ^2^2 _ 4.,/,, 
At 52 = one finds 

A(s', S, = 1 7771, 7771, 7712 ) = 7rPps(s, mi, 7772) <5(s' " s), (34) 

and 

oo 

F'='(0) = (ei+e2) j ^^^^^Pps(s,mi,m2) = ei + e2 (35) 

(mi+m2)^ 

As we have pointed out in the Appendix A, this is just the Ward identity consequence. 

To reveal the relationship between the dispersion integral (^) and the light-cone technique, we introduce 
the light-cone variables 

fc- = -^(fco-fc.); fc+ = -i=(fco + A:.); fc^ = 2fc+fc_ - ^i; (36) 
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into the integral representation for the form factor spectral density (|3l|). We choose the reference frame in which 

- 0, q+ = 0, ql = -q^ 
that is possible at q^ < 0. Performing fc_ integration and setting {fi = +) in both sides of ( p2| ) one finds 

A(s ,s,g |mi,TOi,m2) = — / — « - "j J (37) 

lOTT J xil — x) \ 1 — X X xyl ^ x) ) 

^^U_^_^^^^^±^\ is' + s-2im,^m,r--!-q^) 
\ I — X X x{l — x) J I — X 

Here we denoted x = k2+/P± and fc^ ~ fc2±- 

Substituting @ into Q and performing s and s' integrations, one derives 

H{q^,mi,m2) — J dx(fk±ip{x, k±)'ip{x, k± + xq±)l3{x, k±, q±) (38) 

where the radial light-cone wave function of a pseudoscalar meson is introduced 



^ ^3/2^8(5 -M2)^2;(l- a;) 1 - x a; x(l - 2;) ^ ^ 

^-(mi-m,)^ + fcW(l-^) (3iq^=0) = l 

\/ s — (mi — m^)"^^ s' — (7711 — 7712)2 

The quantity /3 accounts for the contribution of spins. It is different from unity at ^ because both the spin- 
nonflip and spin- flip amplitudes of the interacting quark contribute. The eq.(^) is the normalization condition 
for the soft radial wave function 

' dxSk^\4,(x,kx},S^ = \. (40) 



In terms of this wave function, the pseudoscalar meson axial-vector decay constant /p is represented as 

/"j^ j2; 1 ^ ^2(1 — x) + miX 



fp ~ 9a— f= — J- I dx d k± ^{x,k±) — ■^^^^^=^^^^= (41) 
v27r-^/2 J V s — (7711 — 7772)2 

This expression can be easily deduced by introducing the light -cone variables into the dispersion representation 
(|lO|), making use of ( [Tl] ) and examining the fJ. — + component of the axial current. 

Similarly, introducing the light-cone variables into (|2^ ) yields the following expression for gp^^y 

m\fWc f dxd'^k± , , I X X 
V27r3/2 J ^x{l - x) [m^ + k^Jy/s 

The same expressions for the form factor, pseudoscalar meson electroweak constant, and the two-photon 
decay constant as (^8|)-(^2[) were derived within the light-cone approach in refs 



3 Form factors of meson transitions 



In this section we examine the electroweak transitions of pseudoscalar mesons. First, we derive the dispersion 
representations for transition form factors at <; g and demonstrate them to be equal to those obtained 
within the light-cone calculations. Second, these dispersion representations allow us to perform the analytic 
continuation and derive the form factors of semileptonic decays of pseudoscalar mesons at > q where the 
direct application of the light-cone technique is hampered by the contribution of pair-creation subprocesses. 

^Our Gv{s) is just equal to ho^P) of rcf.[|7| . 



8 




Figure 5: The dispersion graph for the decay < P2|V^^"|^'! >■ 
3.1 The pseudoscalar meson transition form factor at < 

The ampUtude of the weak transition of pseudoscalar mesons Mi M2 (Fig.H) is determined by the two form 
factors F4- and F_ 

< PM,,M2\VnPAh ,Mi> - (Pmi + Pm,)F+{s3) + {Pm^ ~ Pm,)F-{s^) (43) 



a( 

-Pm2 = ^'^h Pi'h ^ ' ^Mi ~ Pm-2 = Pms 7 Pll^ = 

The weak currents are defined through current quarks 

^r = 9i(0)>92(0), = 92(0)7^7591(0). (44) 

The structure of the mesons is described in terms of the constituent quarks by the vertices 

. Q2M^5£3(zMG.,(Pf), M2 : Mi)!^^£i(zMG„2(P|) (45) 

For calculating the tranition amplitude (^3|) we again need the constituent quark matrix element of the weak 
current which is taken in the form 

< Ql(fcl)|qi(0)7,,q2(0)|g2(fc2) 0l(fcl)7pQ2(fc2)/2l(s3) (46) 

The dispersion representation for the form factors reads 

F±(s3) = /2i(.3) y ,(,^ _ Ml) .{S2 - Ml) ^^(^1' ^3|mi, m2, ma) (47) 
Here A± are the double spectral densities of the Feynman graph corresponding to Fig.^ in si~ and S2— channels 

/ dkidk2dhS{kl ~ ml)S{kl - ml)6{kl - ml)5{Pi - k2 - k3)S{P2 - k^ - fci) 
Stt J 

X Sp (^(fci + mi)7^(fc2 + m2)i75(m3 - ^3)175) = (Pi + P2)^A+(si, S2, S3|mi, 7712,7713) (48) 

+ (^^1 - ^2)pA_(si,S2,S3|777i, 7772, 7713) 

with Pi = P2 + P3, Pi = si, Pf = S2, P3 = Pms = '*3- The vertices Gvi^2 are normalized in accordance with 
(H). The trace reads 

Sp (^{ki + 7ni)7p(fc2 + "72)75(7713 - 1^3)75) (49) 

= 2fci^(si - (7772 - 7^,3)^) + 2fc2p(s2 - (7773 - 777i)^) + 2^3^(53 - (mi - 7772)^) 
= (Pi + P2 - 2/03)^ (a(si, m2, 7773) + a(s2, 7773, 7771) - a(s3, 1711,1712)) 
+ (Pi + P2)^a(s3,f«i,"^2) + (^^1 - P2)t, (a(s2,^7^3,"^l) - 0(31,7772,7773)) 
where a(s, ^1, ^2) = s - (/ii - A*2)^- 



<Pm,,M2\AI''\Pm,,Mi> = 
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Making use of the relation 



with 



A(Sl,S2,S3) A(Si,S2,S3j 



6+(si,S2,S3) = -S3{si + S2 - S3 + mj + ml - 2ml) - {ml ~ ml){si ~~ S2) (51) 
6_(si,S2,S3) = {ml - ml){2si + 2s2 ~ S3) - {si - S2){si + S2 - S3 + ml + ml- 2ml), (52) 

we come to the following result for A± 

A±{si,S2,S3\mi,m2,m3) = ^%^^i^^^^^A(si, S2, S3IW1, m2, TO3) (53) 

A(Sl,S2,S3) 

B+(si,S2, S3) = &+(si,S2, S3)(a(si,m2,r7i3) + a(s2,m3,TOi) - a(s3, mi, TO2)) + 0(53, toi, m2)A(si, S2, S3) 
-B_(si, S2, S3) = 6_(si, S2,S3)(a(si,TO2,m3) + a(s2,TO3,mi) - a(s3,TOi,m2)) 
+ (a(s2,m3,TOi) - a(si,m2,m3))A(si,S2,S3) 

3 . 



Here A, the double spectral density in si and S2— channels of the Feynman triangle graph r{Pi , P2 , P3) with 
scalar constituents, is introduced 



y.(p2 p2 p2N _ i f dkidk2dk3d{Pi ~ fca ~ fc3)'^(P2 ^ ^3 - ki) 

^^' 2 '^3 J (2^)4 y (^2_^2_io)(m2-fc2_^0)(m2-fc2_iO) ^ ^ 

= / / ^^\,f2A / 2N ^(^l'^2,S3|mi,m2,TO3). 

7 7r(si - Aff ) 7r(s2 - A/|) 
At S3 < this spectral density reads 

I . 6'(-?'+(si,S2,S3)- A(si,S2,S3)A(s3,mf,TO^)) 

A(si, S2, S3 mi, m2, 7713) = T7TT727 \ ' ^3 < 

16AV^(S1,S2,S3) 

The solution of this 6'-function reads 

S2>(mi+m3)^, sj"(s2, S3) < si < s;i'"(s2, S3); (55) 

X ^S2S3 — S2(mi + m2) — S3(mi + mg) + {ml — m\){ml — m^) ± A^/^(s2, mg, mi)A^/^(s3, mj, m2)^ 
The final dispersion representation for the form factors at S3 < takes the form 

CXJ S+(S2,S3) 

^ I ^ r r \ f ds2Gy2{s2) f dsiGyi{si) B±{si,S2,S3) 

P± S3 =/2i(s3 / TjaT / -7 ,.2^ x3/2r V 56) 

J 7r(s2-M|) J 7r(si - Mf) A3/2(si,s2,S3) 

This representation will be the starting point for the consideration of the meson decays in the next section. 

To demonstrate the equivalence of the dispersion method and the light-cone approach, we turn back to 
the equation (^8|) and again make use of the light-cone variables (|3^), choosing the reference frame PM3+ — 
0, Pmi± = 0. Setting A* = + and making use of gives for A+ 

A+(si, S2, S3|mi, m2, m3) = / — !— ^^d^fc3j_(5(a;i - 0:2)^(1 - xi - X3) (57) 

IbTT J X1X2X3 

( ml ml kl±\^f ml ml {k3^+X3P3±_f 
xd si ] \ S2 



X2 X3 X2X3 J \ Xi X3 X2X3 

(Hereafter Xi — ki+/P+,Pi+ = P2+ = P+,q± = P3±, —q± — S3,x = X3,k± = k3±.) Substituting (57) into ( ^ ) 
yields the following expression for the form factor F+ which gives the main contribution to the semileptonic 
meson decay rate 

dxd^kx^ Gi,i(si) Gi,2(s2) 



F+{q^) - /2i(gJ / i67r3a:(l-x)4si-M2),r(s2-M|) ^^^^ 
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X ^si + S2 - (mi - ma)^ - (m2 - wa)^ + ^ ^ ^ l*^! ^ '^2)^)^ 

Introducing the radial light-cone wave funcion according to ( |39| ) leads to the familiar light-cone expression 

(cf-i) 

F+{q±_) = f2i{q±) / dx(fk±tl;i{x,k±)i{;2{x,k±+xq±)/3+{x,k±,q±); (59) 



_ si + S2 - {mi - 7713)2 - (7772 - 7773)2 + ^(-q2 _ (^^^ _ ^^^2) 
P+ — , = — , = 

2VSl - ('«2 - ^,3)2^52 - (7773 - 777l)2 

_ (7771a; + 7773(1 — a;))(7772a; + 7773(1 - x)) + k±{k± + xq±) 

X{1 - x)^/3l- (7)7.2 - nT.3)2^S2 - (7)7.3 " TOl)^ 

3.2 The transition form factors at > 

For the description of decay processes the form factors in the region < S3 < (Mi — M — 2)^ are necessary, 
while the light-cone representation ( |59| ) is valid only at S3 < 0. For deriving the form factors at S3 > the 
dispersion representation ( |56| ) turns out to be a convenient starting point. We write this representation in the 
following form 

N f ( \ f dsiG^i{si) ds2Gi,2(s2) J5(si,S2,S3) 
i^ S3) = /21(S3 / —, r72T-7 iT^TTT -A(si, S2, S3 7771, 7772, 7773 60) 

J 7r(si - Mf ) 7r(s2 - M|) A(si,S2,S3) 



where A is the double spectral density of the Feynman graph F with scalar constituents (54). This double 
dispersion representation defines the analytic function of S3 both at negative and positive values provided the 
proper expression for the spectral density A is used. It is important to point out that the functions Gt,(s) have 
no singularities in the right-hand side of the complex s-plane , and B and A are polynomials. So the details 
of the dispersion integration at S3 > are determined by the behavior of the quantity A. 

A detailed consideration of the double spectral density A for two massless constituents was performed in 
pOl. We extend that consideration to the case of arbitrary nonzero masses. The same analysis of A for arbitrary 



masses was done by Azimov |30 



Following 0, we first consider the single dispersion relation in A standard calculation yields 

00 

r(p2,P|,p2)= j -^^^^a,{Pls,,Pil (61) 

(mi+m3)2 

where 

a2{Pl S2, P3') = <y+{Pl S2, P|) - <y-iPl S2, (62) 

CT±(S1, S2, S3) = r log (-S2(si + S3 - S2 + 7772 + 7773 - 27772) " i^i ~ s.i){m\ - m^)) . 

167rA(si, S2, S3) 

Hereafter we assume 7772 > 7771. The single dispersion representation reproduces the exact value of the Feynman 
expression (^). Next, we consider the function 0-2(^2, S2, P|) as the analytic function of si = at fixed S2 
and S3 = P| > 0. As S2 < S2 such that 

/ / „_ I ™2 ™2\ 2 

S3 < (7772 - 777i)2, (63) 

both of the functions fT+ and (t_ have square-root branch points on the physical sheet at sf = {^/s2 ~ a/^s)^ 
and sf = {^/s2 + ^/ss)^, connected by the cut (dashed line in FigJ^a). 

In addition, the function has a logarithmic cut on the physical sheet from Si to defined by the 
expression ( |55| ) . The square-root cuts cancel in (T2 = (t+ — ct- , and the logarithmic cut is the only singularity 
of tT2 on the physical sheet. The function (74. has also a logarithmic cut from sj~ to sf which is located on 
the second unphysical sheet of the Riemann surface of the square-root (dotted line in Fig.|^), and does not 
influence the double spectral density. The situation changes at S2 — s" which is determined by the condition 
^i'{^2) — ^r('^2)- The logarithmic and square-root branch points coincide, and for further increasing S2 > S2 
the logarithm branch point moves up through the square-root cut onto the physical sheet, whereas the position 
of the logarithm branch point of (T_ goes to the second sheet (Figj^). Hence, on the physical sheet the function 
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Figure 6: The location of the singularities of (T2 at S3 > 0: a). S2 < b). S2 > Sj. 

(7+ acquires the logarithmic cut from Si to Sj^, and cr_ still has the logarithmic cut from Si to s^. Both of 
the functions have also square-root branch cuts from to s^. In the difference (72 = 0"+ — cr- the square-root 
cuts cancel each other, but the logarithmic cuts add. The resulting expression for the double spectral density 
takes the form 

, I , 0(32 - (mi + 1713^)6(3^ < Si < s+) 

A(si, S2, S3 mi, m2, ms) = T^TTTTf ^ 

16Al/^(Si,S2,S3) 

2g(s3)g(s2-s°)g(sf <si <sr) 

16Al/2(si,S2,S3) 

One can check the double dispersion representation (Q) with the spectral density A given by (^) to reproduce 
correctly the Feynman expression. The first term in (|6^ ) relates to the Landau-type contribution emerging 
when all intermediate particles go on mass shell, while the second term describes the non-Landau contribution. 

In addition to the quantity A, the spectral density of the representation ( |60| ) involves the factor 1/A(si, S2, S3) 
which is singular at the lower limit of the integration in the non-Landau term, namely 

A(si,S2,S3) - (s-sf)(s-sf). 

As it has been discussed in |l0|| , in this case an accurate application of the Cauchy theorem yields the subtracion 
term in the non-Landau contribution. Representing CT2 as a contour integral, we must take into account the 
nonvanishing contribution of the small circle around the point 3i . Underline once more that the presence of the 
factor Gt,i(si) does not change the argumentation as the function G't,(s) has no singularities at si > (m2 +m3)^. 
The final properly regularized representation for the form factors at < S3 < (m2— mi)^ takes the form (omitting 
the constituent transition form factor /2i(s3)) 

00 ^1 

P( \ _ f ds2Gy2{32) f dsiGyijsi) g(si,S2,S3) 

J 7r(s2-M|)y 7r(si-M2)16A(si,S2,S3) ^ ' 



+ 29(33) 



d32Gy2{s2) f d3iGyi{3i) 

7r(s2 - Af|) J 167r(si - sf )3/2 



Gi,i(si)B(si, S2, S3) G't,i(sf )B(sf , S2, S3) 



(si - sf )3/2(si - Ml) (sf - sf )3/2(sf - Ml) 



It should be pointed out, that although the representations (^) and (^5|) were deduced for the case of pseu- 
doscalar mesons, transition form factors of any hadrons have the same structure. A particular choise of the 
initial and final hadrons yields a specific polynomial B. So the performed analysis is valid in the general case 
of hadron decay. 
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4 Calculation results 



We are now in a position to apply the developed formalism to the analysis of the properties of pseudoscalar 
mesons and to the direct calculation of the decay form factors. To this end we must specify the parameters of 
the model, i.e. input the vertex functions of the pseudoscalar mesons and constituent quark masses. 



4.1 Parameters of the model 

For a pseudoscalar meson built up of quarks with the masses mq and m^, it is convenient to introduce the 
function (/) related to the vertex function as 



(m|-m2)2 ^^^2 ^l/2( 2 2) 

The normalization condition (pj) for Gy yields the following normalization condition for 

(f>^{k)k^dk = 1. (67) 



The function ip is the ground-state iS-wave radial wave function of a pseudoscalar meson for which we choose 
a simple exponential form 



4){k) = exp ( -4a— ) (68) 



where fip = rnqniql {mQ + mq) is the reduced mass. The parameterization (|68| ) is inspired by the nonrelativistic 
quantum mechanics and, as we shall see later, is convenient for the analysis of the case mg oo. 

In the nonrelativistic quantum mechanics a bound-state wave function is determined by the motion of the 
particle with the mass /ip in the potential independent of masses, and thus a does not depend on the masses 
as well. Relativistic effects destroy this simple feature of the wave function. In QCD the situation is much 
more complicated because additional dimensional quantities such as Aqcd and the condensates appear. So, a 
should be considered as some unknown function of the quark masses. It is possible to obtain the information 
on the behavior of a as a function of mg at fixed mq = m„^d = 0.25 GeV in the two regions: at small mg and 
niQ — !■ oo. 

At niQ < 0.5 GeV the value of a can be determined by describing the data in the light-meson sector. The 
light-quark masses given in Table ^ and a^^ ~ ctR — 0.02 provide a good description of the data on /jr, /x, and 
the elastic form factors (Figs. |l^ and|l^). The meson decay constants and form factors are calculated with the 
values g\{AP) = 1 and fc{q^) = /c(0), respectively. 

In the region mq — > oo the behavior of a can be found on the basis of the heavy quark symmetry. To this 
end, let us consider the amplitudes of the elastic and inelstic transitions between pseudoscalar mesons consisting 
of heavy Q and light q quarks and introduce the dimensionless form factors as follows 

<M,P'\Qj^Q\M,P>=iP' + P)F,i{q^); < (69) 

F,iiq^) = Kiiuo) ^l~pli(u,-l) + 0{{u; ~ if), u; ^ 1 - ^ > 1 
< M2, P2\Q2lt.Qi\Mi,Pi >= [Pi + P2)F+{q^) + (Pi - P2)F_{q^)- < q^ < (Mi - Ahf (70) 

^,( -) -^;^^-,(:^) , ^^^-(:-): . = M^±M^>i. 
^ ' 2VM1M2 ^ ' 2VM1M2 ^ ' 2M1M2 

h+{uj) = h+{l) - p''{lu - 1) + Oiiuj - if). 

In the limit of infinitely heavy quarks Qi,2, the amplitudes are expressed in terms of the single universal 
Isgur-Wise function (IW) ^{uj) [|| 

h+{io)^hM^aco), h^{uj)^0, ^{uj)^1-p\lu-1) + 0{{u;-1)^). (71) 

In addition, the qeavy quark symmetry predicts the universal relation for heavy-meson decay constants 

^Mp/g = const. (72) 
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within the HQET El 



The asymptotic relations ( |7l| ) and ( |72| ) are the zero-order terms of the l/mg-expansion which is calculable 



A particular form of the IW function depends on the heavy meson wave function. 



The expressions (fTj) and ( [72| ) mean that the HQ symmetry restricts the possible behavior of the meson 
wave function at large mq. Table ^ gives the results on fp and p\i vs mq at ruq — 0.25 GeV ^ and Fig.|| presents 
the quantity ^JrrTqfp as the function of mq for various values of a. In the HQ limit, for a finite binding energy 
of the meson the heavy meson and the heavy quark masses coincide, Mq/mq = 1. So, the value of ^Jmqjp 
should be independent of the heavy quark mass. 

These results show that the asymptotic relations ( [7l| ) and (^2|) are satisfied if the parameter a of the wave 
function ( |68[ ) tends to a constant ocao as mq oo. 

Thus, the function a{mq) has the following behavior: it is equal to 0.02 at mq < 0.5 GeV and tends to a 
constant a^o as mq — > oo. For investigating the B and D mesons and their decays we need the information on 
a in the region mq = 2-^5 GeV. 

The simplest way is to extract a at mq = 2-^5 GeV from the analysis of fpi and fs as we have done 
for the light mesons. In the absence of the experimental data we refer to the results of other models. As one 
can see, the decay constants fp calculated with a from the range 0.02 < an, cub < 0.04 cover the regions 
160 MeV < fn < 230 MeV and 130 MeV < fs < 200 MeV which include the predictions of most of the 
models. Hence, the values of ao and as related to the true wave functions of D and B mesons are expected to 
be inside the interval 0.02 0.04. 

However, there is an attractive possibility to specify ajj^s more precisely. Namely, it seems reasonable to 
assume a to be approximately constant in the region mq > l-;-2 GeV . There are at least two arguments behind 
this assumption. Firstly, a system consisting of a heavy and a light particles behaves like a quasinonrelativistic 
system. And secondly, there are no visible sources within QCD to yield steep changes of a in this region. Then 
for the B and D mesons one expects ao = olb — Ooo- The next step is to estimate aoo- We consider the 
value Qfoo — 0.02 to be both attractive and reasonable: on the one hand, the same parameter describes all 
ground-state mesons, and on the other hand, one finds for ctoo = 0.02 



JPao 



5.8 Gel/3/2 



in agreement with the value 0.6 0.7 estimated in [Q. 

Assuming ao = cub — ctoo, we can estimate the magnitude of the higher order 1/mq corrections which 
determine the deviations of the calculated fp and p^; at finite mq from the asymptotic relations (|7^) and (|7^) . 
Rather strong violation of the HQ symmetry for b— and c— quarks (5 -f- 15% at mq = 5 GeV and 20 30% at 
mq = 2 GeV) can be observed both in fp and p^i at aoo = 0.02 0.04. 

We shall analyze the transition form factors obtained at ajj^B = 0.02 and 0.04. If our assumption an = 
ub — ckoo = 0.02 does not work properly, the form factor calculations for a = 0.02 and a = 0.04 give an interval 
which is expected to include the true value. 

Table | gives the numerical parameters of the model. 

4.2 Discussion 

1. The results on the axial-vector decay constant fp are shown in Fig.^ and Table ^. Assuming a{mq) = aoo at 
mq > 2 GeV, one can see the asymptotic relation ^ymqfp = const to work perfectly at mq > 40 — 50 GeV, and 
finds essential corrections to the asymptotic relations at lower mq. For aoo = 0.02 one obtains fo — 234 Mel^ 
and /b = 202 MeV that confirms the expectation fn — fB @. These values for the decay constants correspond 
to the constituent quark decay constant 5^ = 1. In reality, the latter can be less than unity, ~ 0.75 1. 
This will lead to decreasing the fp. 

2. Figures 12-17] present the elastic and transition form factors calculated with a^.B ~ 0.02 and 0.04. 



The if — > TT transition form factor is well approximate by the linear function F^{q'^) = -F+(0) + aq^, 
F+{0) = 0.96, a = 1.27 GeV'^ in agreement with the results of (l). 

The parameters of the monopolc F+iq"^) = F+(0)/(l - q'^/M^^J and the dipole F+{q'^) = F+(0)/(l - 
l^dipY fits to the other transition form factors are given in Table The dipole formula excellently ap- 
proximates the transition form factors with better than 1% accuracy. Although the monopole fit provides a 
worse accuracy, its parameters agree with the vector meson dominance. The values F+(0) are close to the 
corresponding results of QCD sum rules (cf. Tables ^ and |^) and the existing experimental data. 
3. Fig.|l^ plots the IW function ^(w) = ft.+ (cj) for the decay B D aX various values of au and as. Table || 
gives the parameters of the calculated IW function. The function [uj] turned out to be negligibly small in 
agreement with ([7l|). 

The IW function has been extensively studied both theoretically and experimentally (see Table 0). The 
analysis by ARGUS js^ and most of the earlier theoretical results suggested 1 < < 2. However, a recent 
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analysis by CLEO as well as recent theoretical estimates favor the lower values < 1. We found the relation 

0. 7 < < 0.9 for all values of aD,B from the considered interval. 

As it follows from the HQ symmetry, the value ^(1) strongly depends on the relationship between and 
as- it turns out to be close to unity for as ^ cub and steeply decreases as as 7^ cud- One can find rather 
uncertain constraint for the considered region of the parameters aD,B 0.87 < ^ 0.98. 

Let us underline that except for the relationship between ao and as, the value ^(1) is also affected by the 
particular values of heavy-meson binding energies. At large quark masses and the binding energy kept finite, 
the heavy meson and heavy quark masses coincide, Moo /moo = 1- Hence, the positions of the 'quark zero recoil 
point' gQ = (mQj^ — mg^)^ and the meson zero recoil point q^^^^ = (Mi — M2)^ also coincide. For infinitely 
heavy quarks this yields (^(1) = 1. For the physical heavy quarks and mesons, the positions of the 'quark zero 
recoil point' ujq — 1 ~ ((?71qi — mQ2f' ~ (^^i — M^Y) /2M1M2 and the 'meson zero recoil point' a; = 1 do not 
coincide any longer. The calculated ^(t^-'o) turns out to be not far from unity if au = as- So, the value ^(1) 
is sensitive to the particular values of the quark masses. The quark masses used in our calculation are chosen 
such that rrih — rric — Mb — Mu, and thus g§ — q^ax ^^'^ ^0 = 1- That is why ^(1) ~ 1 at an — olb- For other 
reasonable values of quark masses, the deviation from unity at an — cub are found at the level of 3-4%. 
4. The analysis of the analytic properties of the hadron transition form factors yields the following typical picture 
demonstrated in Fig.|6[ at < the contribution of the non-Landau singularity is absent, and the Landau- 
type singularity determines the form factor; in the region < < (to2 — mi)^ both of them are essential; at 
the point — (m2 — mif' the contribution of the Landau singularity vanishes, and the non-Landau singularity 
determines the decay form factor at this 'quark zero recoil' point. 

For hadron decays related to the heavy-to-heavy quark transitions, a specific relationship between the 
Landau and the non-Landau contributions to the dispersion representation is observed: the normal Landau 
contribution dominates the form factor at all q^ < (to2— mi)^, whereas the anomalous singularity is essential only 
in the close vicinity of this point. So, effectively the transition form factor are determined by the contribution 
of the Landau contribution only. Thus, the HQ symmetry can be formulated in the language of the analytic 
properties of the transition form factors as the dominance of the Landau singularity in the almost whole 
kinematical region. 

In the case of the meson decay related to a heavy-to-light quark transition, the anomalous non-Landau 
contribution is important in a broad kinematical region. So the relations suggested by the HQ symmetry would 
not work properly. 

5 Conclusion 

We investigated form factors of hadron transitions within the relativistic constituent quark model and proposed 
a formalism for a direct calculation of hadron decay form factors. The developed approach was applied to the 
analysis of the electroweak properties and transitions of pseudoscalar mesons. Our main results are: 

1. The equivalence of the light-cone constituent quark model and the approach based on the dispersion relation 
integration over a bound state mass for the description of leptonic decays and transition form factors at spacelike 
momentum transfers has been demonstrated. Although the comparison has been performed for a particular 
case of pseudoscalar mesons, the approaches are equivalent for the description of any hadrons. 

2. The obtained dispersion formulation of the light-cone constituent quark model allows a consideration of the 
decay processes where the direct application of the light-cone technique is hampered by the contribution of 
pair-creation subprocesses. The analytic continuation in the dispersion representation of the transition form 
factor yields the form factor at timelike momentum transfers expressed through the meson radial light-cone 
wave function. Along with the normal Landau singularities, the anomalous non-Landau singularities contribute 
to the form factors at q^ > 0. 

3. For hadron decays related to the heavy-to-heavy quark transitions a specific relationship between the 
contributions of the Landau-type and the non-Landau singularities has been observed. This allows a formulation 
of the heavy quark symmetry in the language of the analytic properties of the decay form factors as the 
dominance of the normal Landau contribution in the almost whole kinematic region of momentum transfers. 

4. Electroweak properties and form factors of pseudoscalar mesons have been analyzed using a parameterization 
of the meson wave function based on the heavy quark symmetry. We have examined the dependence of the 
axial-vector decay constant on the heavy-quark mass, and found fo — 235 MeV and fs — 200 MeV. These 
values can be decreased by a factor of 0.75-^ 1, if the decay constant at the level of the constituen quarks is less 
than unity. 

The correlation between the axial-vector decay constant fp and the transition form factors yields the IW 
function parameter p^ = 0.8 ± 0.1 for the axial-vector decay constants from the intervals 160 MeV < /_d < 
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235 MeV and 130 MeV < fo < 200 MeV 

Analyzing the dependence of fp and the heavy meson form factor on the heavy quark mass we have found 
that the violation of the HQ symmetry relations can be expected at the 10-20% level for the b- and c-quark 
masses . 

5. The calculated form factors of pseudoscalar meson transitions have been approximated with a 1%-accuracy 
by the dipole formula in the whole kinematic region. The form factors are also compatible with the vector 
meson dominance and are close to the results of the QCD sum rules. 

The developed approach can be applied to the description of the pseudoscalar-to-vector meson transitions 
and rare decays of heavy mesons. This work is now in progress. 

I am grateful to V.V.Anisovich, Ya.I.Azimov, and K.A.Ter-Martirosyan for discussing the general problems 
and technical details related to hadron decays. I am also indebted to the German Ministry of Science and 
Technology for the financial support at the early stage of this work and to H.R.Petry for his hospitality during 
my stay in Bonn. 
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Figure 7: One of the terms in the expansion of Ao{s) 



6 Appendix A: Bound state description within dispersion relations 

To ihustrate main points of the dispersion approach we consider the case of two spinless constituents with the 
masses mi and m2 interacting via exchanges of a meson with the mass fi. We start with the scattering amphtude 

A{s,t)=<k[,k!,\S\ki,k2>, s^{ki + k2f, t={ki-k[f (73) 

The amphtude as a function of s has the threshold singularities in the complex s-plane connected with elastic 
rescatterings of the constituents and production of new mesons at 

s = {nil + 7712)^, (mi + m2 + fi)^ , (mi + m2 + 2//)^ . . . (74) 

We assume that an 5- wave bound state with the mass M < mi +m2 exists, then the partial amplitude Ao{s) has 
a pole at s = Af ^. The amplitude A{s, t) has also i-channel singularities at i = {njiY; n ~ 1,2,3... connected 
with meson exchanges. If one needs to construct the amplitude in the low-energy region s > (mi + m2)^ the 
dispersion N / D representation turns out to be convenient. Consider the 5- wave partial amplitude 

1 

Ao{s) J dzA{s,t{s,z)), (75) 
-1 

where t{z) = —(1 — z)X(s,m1,m2)/2s, z — cos6 in the c.m.s. The Ao{s) as a function of complex s has the 
right-hand singularities related to s-channel singularities of A{s,t). In addition, it has left-hand singularities 
located at s = (mi + m2 )2 - (nn)^; n = 1,2,3.... They come from i-channel singularities of A{s,t). The 
unitarity condition in the region s w (mi + m2)^ reads 

Im^o(^) = p(^)|^o(.)^ (7g) 

with p{s) the two-particle phase space. The N/D method represents the partial amplitude as Ao{s) = 
N{s)/ D{s), where the function N has only left-hand singularities and D has only right-hand ones. The unitarity 
condition yields 

00 

f dS p(S)N(S) 

D{s) = l- / — ^ ^ = 1 ~ B{s). (77) 

J TT S — S 

(mi+m2)^ 



Assuming the function N to be positive we introduce G{s) — ^ N{s). Then the partial amplitude takes the 
form 

Aois) = G{s) [1 + B{s) + B'{s) + B'is) + . . .] G{s) = (78) 
This expression can be interpreted as a series of loop diagrams of FigJ^ with the basic loop diagram 

Bis)= J ^m^. (79) 

J TT S — S 

The bound state with the mass M relates to a pole both in the total and partial amplitudes at s = so 
B{M'^) = 1. Near the pole one has for the total amplitude 

A = < k[,k2\P > < P\ki,k2 > +regular terms 

^ XUklk'2) ^2^ p2 XP{kl,k2) + ■ ■ ■ (80) 
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Figure 8: One of the terms in the series for T^. 

where Xp(fci, fc2) is the amputated Bethe-Salpeter ampUtude of the bound state. The dispersion amplitude near 
the pole reads 

A — N/ D + regular terms related to other partial waves 



(M2 - s)B'{]VP) 



(81) 



where Gy is a vertex of the bound state transition to the constituents. The singular terms correspond to each 
other and hence ^ 

Xp(fci, fc2) -> G„(P') = \^ (82) 

Underline that among right-hand singularities the constructed dispersion amplitude takes into account only the 
two-particle cut. 

Let us turn to the interaction of the two-constituent system with an external electromagnetic field. The 
amphtude of this process =< k[, kl2\Jf_t{q)\ki, k2 > in the case of a bound state takes the form 



T, = <k[,k',\P'>-——-<P'\J,{q)\P>^^^-^,<P\klM> 



XpiK.k'^) p,^^ ^^ {P' + P)^F{q^)p^^j^Xp{kiM) + . . . (83) 



where the bound state form factor is defined as 

<P'\J^{q)\P>={P' + P)^F{q^) (84) 
The dispersion amplitude with only two-particle singularities in the P^- and P'^-channels taken into 



account is given |28| by the series of graphs in Figj8| 

These graphs are obtained from the dispersion scattering amplitude series by inserting a photon line into 
constituent lines. The amplitude reads 

T^(P', P, q) = 2P^{q)T{s', S, q') + ^C, (85) 

p2 = s, P'2 = s', q = P'-P, PM = - ^ '?)^' 

The dispersion method allows one to determine r(s, s', q^), which is the part of the amplitude transverse with 
respect to q^. Summing up the series of dispersion graphs in Fig.H gives 

Here 

r(s ,s,9 ) = / —7Z r—-r, rA(s , S, g ), 



7r(s — s) 7r(s' — s) 

and A(s', s, g^) is the double spectral density of the three-point Feynman graph with a pointlike vertex of the 
constituent interaction. 

The longitudinal part C is given by the Ward identity 

^-T^W')-B(.))^ (87) 
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At s = s' = M^, the quantity develops both s and s' poles, so 



T^iP', P, q) = _ j {P' + P)^F(g') j + less singular terms 



ne)= I 4^^4^;^^A(.',.,,^). (89) 



where 

oo 

P(n-2\ - [ 

7r(s - KP) 7r(s' - A'P) 

is the bound-state form factor (see ( ^2|) and (|8^)). So, the quantity < P'| J^((?)|P > corresponds to the three- 
point dispersion graph with the vertices G„. The following relation is valid A(s',s,0) — tt6{s' — s)p{s). This 
is a consequence of the Ward identity which relates the three-point graph at zero momentum transfer to the 
loop graph. This relation yields the charge normalization F{0) = 1. The expression ( ^9|) gives the form factor 
in terms of the A^-function of the constituent scattering amplitude and double spectral density of the Feynman 
graph. In general, the following prescription works: to obtain the dispersion expression spectral density in 
channels corresponding to a bound state, one should calculate the related Feynman graph spectral density and 
multiply it by Gy. 

If the constituent is a nonpoint particle, the expression ( |89| ) should be multiplied by form factor of an on-shell 
constituent. 
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Figure 9: The quantity niQ^fp as the function of mg at — 0.25 GeV. 



Table 1: The form factors of the decays B tt, p a.t — 0. The labels Lat, SR, QM and LCQM stand for 
Lattice, Sum Rules, Quark Model and Light Cone Quark Model, respectively. 







F+(0) 


V{0) 


Ai(0) 


^2(0) 


y(o)Mi(o) 


v42(0)Mi(0) 


Lat 


1 
1 
1 
I 


16 
16 
17 
17 


a 
b 
a 
b 


0.29± 0.06 
0.35± 0.08 
0.26± 0.16 
0.30± 0.19 


0.45±0.22 
0.53±0.31 
0.34±0.10 
0.37±0.11 


0.29±0.16 
0.24±0.12 
0.25±0.06 
0.22±0.05 


0.24± 0.56 
0.27± 0.80 
0.38± 0.22 
0.49± 0.26 


2.0±0.9 
2.6±1.9 
1.4±0.2 
1.6±0.3 


0.8±1.5 
1.0±3.1 
1.5±0.7 
2.3±0.9 


SR 




10 
12 




0.24± 0.025 
0.40± 0.20 












QM 


WSB|gl 
GlSWg 


0.33 
0.09 


0.33 
0.27 


0.28 
0.05 


0.28 
0.02 


1.2 
5.4 


1.0 
0.4 


Table 2: The form factors of the decays D K, K* at ^ 0. 






F+(0) 


^(0) 


^i(O) 


A2{0) 


^(O)Mi(O) 


^2(o)Mi(o) 


Exp 


[ 


24 




0.77± 0.04 


1.16±0.16 


0.61±0.05 


0.45± 0.09 


1.90±0.25 


0.74±0.15 


Lat 


[ 
[ 


16 
17 




0.78± 0.08 
0.60± 0.22 


1.08±0.22 
0.86±0.24 


0.67±0.11 
0.64±0.16 


0.49± 0.34 
0.40± 0.32 


1.6±0.3 
1.3±0.2 


0.7±0.4 
0.6±0.3 


SR 


[ 


11 




0.6± 0.15 


1.1± 0.25 


0.5± 0.15 


0.6± 0.1 


2.2±0.2 


1.2± 0.2 


QM 


WSBpl 
GISW[| 


0.76 
0.8 


1.23 
1.10 


0.88 
0.80 


1.15 
0.80 


1.4 
1.4 


1.3 
1.0 


LCQM 


li 


0.73 


0.92 


0.63 


0.42 


1.46 


0.67 
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Table 3: The decay constants fp of pseudoscalar mesons MeV. 





TT 


K 


D 


B 


Exp |2| 


130.7± 0.46 


159.8±1.9 


< 310 




Lattice 19| 






200± 30 


180± 40 


Sum Rules 






160 [11 
165^1951 L4 


130-^ 200jl^ 


LCQM § 
LCQMg 
This work 


130.7 
130.7 
130 


162 
162 
160 


220 
206 
234 


188 
186 
202 



Table 4: The decay constants fp of pseudoscalar mesons built up of quarks with the masses mg and and 
the slope of hei at = 1 calculated from < MqlQ'^f^QlMq > as functions of mq at niq = 0.25 GeV. 





a = 0.01 


a = 0.02 


a = 0.04 


a = 0.08 


mq, GeV 


fp,MeV 




fp, MeV 




fp,MeV 




fp,MeV 




0.25 


151 


0.04 


130 


0.06 


104 


0.08 


80 


0.1 


0.4 


190 


0.25 


160 


0.35 


128 


0.5 


97 


0.65 


1.8 


324 


0.6 


234 


0.65 


163 


0.82 


110 


1.0 


5.2 


308 


0.75 


202 


1.0 


132 


1.05 


85 


1.1 


10 


254 


1.0 


162 


1.05 


102 


1.1 


64 


1.25 


20 


195 


1.0 


122 


1.1 


76 


1.23 


48 


1.45 


40 


143 


1.0 


89 


1.11 


55 


1.25 


34 


1.66 


80 


103 


1.0 


63 


1.11 


39 


1.25 


24 


1.66 



Table 5: The constituent quark masses and the calculated fp for a = 0.02. 



quark 


quark mass,Gey 


meson 


meson masSjGeT^ 


fp.MeV 


u,d 


0.25 


7r+(w(i) 


0.14 


130 


s 


0.40 


K+{us) 


0.49 


160 


c 


1.80 


D+{cd) 


1.87 


234 


b 


5.20 


B+{uh) 


5.27 


202 



Table 6: The parameters of the monopole and dipole fits to the form factor. The masses of the lowest 
vector mesons which are expected to dominate the form factors are given in brackets. 





CtB = OtD = 0.02 


as ~ ctD ~ 0.04 


Decay 


F+{Q) 


Mmon, GeV 


Mdip, GeV 


F+(0) 


Afrf.p, GeV 


B ^ D 


0.73 


5.7 


7.7 


0.68 


7.20 


B ^TT 


0.23 


5.2 [5.324] 


6.2 


0.22 


6.08 


D ^ K 


0.70 


2.22 [2.11] 


3.0 


0.70 


2.95 


D^TT 


0.55 


2.1 [2.01] 


2.8 


0.59 


2.68 



Table 7: The slope of_the IW function p\ 





ARGUS |3ll 


CLEO [|2| 


Lat [|l§ 


SR[1§ 


srH 


This work 




1.07±0.17 


0.87 ±0.12 


1.2 


> 1.04 


0.7 ±0.2 


0.8 ±0.1 



Table 8: The parameters of the IW function. 





as 


^(1) 


P' 


0.02 


0.02 


0.98 


0.78 


0.02 


0.04 


0.87 


0.75 


0.04 


0.02 


0.93 


0.7 


0.04 


0.04 


0.98 


0.88 
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Figure 10: The 7r+ form factor, a,r = 0.02. 
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Figure 11: The K+ form factor, = 0.02. 
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Figure 12: The D+ form factor. Solid - ud = 0.02, dotted - an = 0.04. 
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Figure 13: The B+ form factor. Solid - as = 0.02, dotted - as = 0.04. 
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Figure 14: The form factor F+{q'^) for B ^ D. Solid - an = as = 0.02, dotted - the monopole fit, dash-dotted 
- the dipole fit. Dashed - an = as = 0.04 




Figure 15: The form factor F+{q'^) for B ^ n. Solid - as = 0.02, dotted - the monopole fit, dash-dotted - the 
dipole fit. Dashed - as = 0.04. 
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Figure 16: The form factor F+ for D ^ K. Solid - an = 0.02, dotted - the monopole fit, dash-dotted - 
the dipole fit. Long-dashed - the Landau singularity contribution, long-dash-dotted - the non-Landau term. 
Dashed - an = 0.04. 
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Figure 17: The form factor for D n. Solid - an = 0.02, dotted - the monopole fit, dash-dotted - the 
dipole fit. Dashed - an = 0.04. 




Figure 18: The Isgur-Wise function ^{w) = h+{co) of the decay B ^ D: SoHd - ao = as = 0.02, dashed - 
Q!£) = as = 0.04, dash-dotted - an = 0.04, as = 0.02, dotted - an = 0.02, ub = 0.04. 
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